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ABSTRACT 

FeLoBALs are a rare class of quasar outflows with low-ionization broad absorp¬ 
tion lines (BALs), large column densities, and potentially large kinetic energies that 
might be important for ‘feedback’ to galaxy evolution. In order to probe the physical 
properties of these outflows, we conducted a multiple-epoch, absorption line variability 
study of 12 FeLoBAL quasars spanning a redshift range of 0.7 ^ ^ 1.9 over rest 

frame time-scales of ~10 d to 7.6 yr. We detect absorption line variability with ^ 8cr 
confidence in 3 out of the 12 sources in our sample over time-scales of ^0.6 to 7.6 
yr. Variable wavelength intervals are associated with ground and excited state Fen 
multiplets, the Mgn AA2796, 2803 doublet, Mgi A2852, and excited state Nin mul- 
tiplets. The observed variability along with evidence of saturation in the absorption 
lines favors transverse motions of gas across the line of sight (LOS) as the preferred 
scenario, and allows us to constrain the outflow distance from the supermassive black 
hole (SMBH) to be < 69, 7, and 60 pc for our three variable sources. In combination 
with other studies, these results suggest that the outflowing gas in FeLoBAL quasars 
resides on a range of scales and includes matter within tens of parsecs of the central 
source. 
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1 INTRODUCTION 

Broad absorption lines (BALs) in rest frame UV quasar 
spectra are typically seen as blue-shifted troughs that indi¬ 
cate the existence of gas outflows. The ions that produce the 
absorption are often highly ionized species such as C IV, N v, 
and O VI; their presence defines the high-ionization BAL (Hi- 
BAL) quasar cl ass, which comprise s between ~10 and 26 per¬ 
cent of quasars dTrump et al.|[200^ ). Sources with additional 
BALs from lower ionized species (e.g. Mgn, Alii, and Aim) 
are categorized as low-ionization BALs (LoBALs) and com¬ 
pose a much smaller fraction of quasars than HiBALs. 

Throughout this paper, we define iron low-ionization 
BALs, or FeLoBALs, as any quasars with Fell absorption 
along wit h transitio ns that characterize HiBALs and LoB¬ 
ALs fsee lHall et al.ll2002al lbl and reference s there in for gen¬ 
eral FeLoBAL information!. iTrump et ahl ([2006) estimated 


FeLoBALs to comprise < 0.3 per cenlQ of quasars. The rar¬ 
ity of FeLoBAL outflows is not well understood; they may 
represent a short-lived phase of quasar evolution, or perhaps 
the Fell-Fem gas covers a very small fraction of the active 
galactic n ucleus (AGN) and only few sight lines permit t heir 
detection (lHamann fe Sabrall2004l ; iHamann et al .1 [20 1 2l ) . 

There are conflicting ideas about where FeLoBAL 
outflows are located, and this has important implica¬ 
tions for the calculated kinetic energies and the possi¬ 
ble role of these outflows in ‘feedback’ processes with 
the host galaxy. FeLoBALs at kiloparsec distances from 
the supermassive black hole (SMBH) may inject enough 

1 The absorbing gas and dust likely produce selection effects in a 
flux-limited sample, and 0.3 per cent may therefore not represent 
the true FeLoBAL quasar fraction. 
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energy and momentum into the int erst ellar medium 
(ISM) to d i srupt s tar formation (|Fanjah__et_al J |2012|; 
iLeighlv et al.l 120141). iFaucher-Giguere fc C'uat aertl (|2012h 
and Faucher-Gieuere. Quataert &: Murray ( 20121 ) provide a 
theoretical framework for the formation of Fell and Fern 
at such radial scales that consists of gas driven from the 
acc retion disk interact i ng with ISM clumps. Oth er mod¬ 
els llMurrav et al.lll995l : Ide Kool fe Begelmanlll995l ) instead 
suggest that all ionic species arise from gas accelerated off 
the accretion disk at sub-parsec scales from the SMBH. 

Studies that have utilized unblended, excited state Fen 
lines and photoionization models have constrained FeLoBAL 
outflows to exist > 700 pc from the SMBH in a handful of 
sources, with estimated kinetic luminosities and mass flow 
rates on the order of ~ l per cent of Lboi and ~100 M^ 


yr 1 , respectively (de Kool et al. 

2001: Korista et al.l 

200|; 

Moe et al. 2009; Dunn et al. 2010 

: Bautista et al. 2010 

). Ki- 


netic luminosities and mass flow rates of these magnitudes 
make these sources viable contributors to feedback and may 
play a r ole in explaining SMBH-bulge relations (e.g. M -a 
relation, iFerrarese fe Merrit3l200d : iGebhardt et al.l[200d ). 

Observational constraints on outflow location and other 
physical quantities can also be extracted using absorption 
line variability studies. Two scenarios that can potentially 
explain variations are that the absorbing gas undergoes an 
ionization change, or there is a change in covering frac¬ 
tion along our line of sight (LOS) due to transverse mo- 
ti ons of gas across the em ission region of the AGN (see §5.1 
of ICapellupo et all 120131 for other possibilities and reasons 
why they are unlikely for most sources). Limits on outflow 
location depend on which interpretation is utilized. An ion¬ 
ization change relates the measured variability time-scale to 
the recombination time of the gas, allowing for an estimate 
of the density, and hence distance from the central source. 
Models invoking absorbers crossing the LOS constrain an 
outflow crossing speed from the measured time-scale, which 
can be related to a keplerian velocity and hence distance 
from the central black hole. 

Work focused on BAL variability associated with high- 
ionizati on species has e stabli shed a range of outflow dis¬ 
tances. iMoe et all (120091 ) and ICapellupo et al.l (120131 ) used 
models involving absorbers crossing the LOS to constrain 
the outflows to be at most sub-p arsec and par sec scales 
from the SMBH, respectively, while iBarlowl d 19941 ) assumed 
an ionization change scenario and established an u pper 
limit of a few hundred parsecs. ICapellupo et al.l (j2013i ) dis- 
cussed the r elative merits of t he tw o interpretations, and 
ICapellupo. Hamann fc Barlow! d2014l ) provided strong evi¬ 
dence supporting a LOS-crossing scenario in the P v BAL 
quasar Q1413+1143. 

Very few published studies of FeLoBAL variability ex¬ 
ist. The properties of the entire FeLoBAL population are 
therefore not well understood due to few constraints on 
the outflows and a lack of statistic a l resu lts from large 
sample variability studies. lHall et all (2011) observed dra¬ 
matic changes in Mgn and a disappearance of the Fen 
lines at the same velocity in FBQS J1408+3054, and con¬ 
strained the outflow to be between 1.7 and 14 pc from 
the SMBH based upon a m odel with clouds moving across 
the LOS. IVivek et al.1 { 20121 ) probed the time variability of 
five FeLoBAL sources, reporting changes from Fern and 
Aim fine structure absorption lines in one of those ob¬ 


jects and concluded that those ions likely originated in the 
ISM based upon ioniz ation considerations and a model from 
iRogerson et al.l (|201ll ). 

We have acquired multiple-epoch spectroscopy for a 
sample of twelve FeLoBAL quasars and report results on 
absorption line variability using objective statistical crite¬ 
ria. Due to the limited work on FeLoBAL variability to 
date, this study provides a significant increase in the in¬ 
formation available on these outflows. Section 2 (§2) details 
the observations and data reduction procedures, §3 explains 
the analytical approach used for investigating variability, §4 
presents our results, §5 discusses the implications of the vari¬ 
ability, §6 summarizes our findings, and §A provides details 
on individual objects of interest. 


2 DATA 


Our FeLoBAL quasar sample repr esents a subse t of objects 
taken from catalogs compiled by iTrumpet al.l tl2006l ) and 
ICibson et al.1 (120091) along w ith J12144 !q, pre viously stud¬ 
ied b y Ide Kool e t al.1 d2002al) . and J142703 (see iBecker et all 
120001) . Trump et al. (120061) identified FeLoBALs by visually 
inspecting their sample of BAL quasars in a redshift range 
0.5 ^ z ^ 4.38, based o n t h e Sloa n Digital Sky Survey 
(SDSS) data release 3. I Gibson et al.l (120091) constructed a 
sample of LoBALs with redshift 0.48 ^ 2 ^ 2.28 from SDSS 
data release 5, and we further identified FeLoBALs from 
visual inspection of those sources. The selection process in¬ 
volved identifying quasars with strong, low-ionization Fell 
or Fem absorption lines and proceeded without any prior 
knowledge of whether a given source was variable. Due to 
selection effects, the conclusions of this study may not apply 
to the entire FeLoBAL population. 

The twelve FeLoBAL quasars in our sample are listed 
in Table 1 with relevant information including redshift, ab¬ 
solute i-band magnitude, dates of observation, wavelength 
coverage, variability, and radio loudness. The range of red¬ 
shift across our sample is 0.7 ^ z ^ 1.9 with time-scales 
between any two observations ranging from ~10 d to 7.6 yr 
in the quasar’s rest framt@. Spectra were taken by the 2.4- 
m Hiltner telescope at MDM Observatory and the 9.2-m 
Hobby-Eberly Telescope (HET) at McDonald Observatory, 
and were supplemented by spectra from SDSS data release 
7 (DR7) and 12 (DR12). We also utilized I/-band magni¬ 
tudes, estimated on the basis of an unfilte red CCD, from 
the Catalina Sky Survey (CSS: |Prakell200'9l ) data release 2 
in our interpretations of the variability (see §5.1.1). 

MDM spectra from 2010 and 2013 were obtained using 
the Boiler & Chivens CCD Spectrograph (CCDSj3 while 
MDM epochs from 2011-2012 employed the Ohio State 
Multi-Object Spectrograph (OSMOS: iMartini et al.ll201lh . 
Slit widths used during observations with CCDS and OS¬ 
MOS were 1 and 1.2 arcsec, respectively, and the slit angle 
was rotated to coincide with the parallactic angle to mini¬ 
mize effects of atmospheric dispersion. Integrated exposure 


2 We will identify sources using the first part of their full names, 
listed in Table 1. 

3 All subsequent time-scales are in the rest frame of the quasar. 

4 http://www.astronomy.ohio-state.edu/MDM/CCDS/ 
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Table 1. FeLoBAL quasar sample. 


Name 

2 o, 

Mi 

HET^ 

DR7 

DR12 

MDM 

Coverage (AW 

Variable?^ 

rJ_ 

J030000.57+004828.0 

0.9003 

-27.6 

— 

2000.73, 2000.92, 
2001.81, (2001.14) 

— 

2012.07 

2000-3130 

N 

0.0 

J031856.62-060037.7 

1.9274 

-28.5 

— 

2001.04 

— 

2011.02, 2012.08 

1300-2030 

N 

— 

J084044.41+363327.8 

1.2253 

-28.4 

2000.10, 2000.20, 
2000.28, 2000.34, 
(2000.23) 

2002.13 


2010.21, 2011.03, 
2012.08 

1930-3260 

T 

4.6 

,1104459.61+365605.2 

0.7010 

-26.8 

— 

2005.25 

2011.15 

2011.02, 2012.07 

2230-5410 

T 

22 

J112220.77+153927.9 

1.1083 

-26.1 

— 

2005.03 

2012.16 

2012.09 

1800-4360 

T 

0.0 

J112526.12+002901.3 

0.8635 

-25.8 

— 

2000.19 

2011.04 

2012.07 

2040-4940 

T 

0.0 

J115436.60+030006.3 

1.3904 

-27.1 

— 

2001.38 

2011.31 

2012.08 

1590-3850 

N 

0.0 

J120337.92+153006.7 

1.2176 

-26.9 

— 

2005.27 

— 

2012.09 

1800-3090 

Y 

5.4 

J121442.30+280329.0 

0.6949 

-26.3 

2000.17, 2000.24, 
2000.33, 2000.46, 
(2000.28) 

2006.24 

2013.18 

2010.21, 2011.31, 
2012.08 

2540-4280 

Y 

14 

J123103.71+392903.7 

1.0038 

-25.7 

— 

2005.26 

2011.08 

2012.08, 2013.05 

1900-4590 

Y 

0.0 

J142703.60+270940.0 

1.1661 

-26.3 

2000.17, 2000.24, 
2000.33, 2000.53, 
(2000.26) 

2006.38 

2012.38 

2010.14, 2010.36, 
2011.31 

1750-4250 

T 

24 

J152438.80+415543.0 

1.2301 

-26.3 

— 

2005.39 

2012.17 

2012.09 

1700-4130 

N 

22 


a Emission line redshifts taken from iHewett fe Wildl feoid ). 

^ Absolute i— band magnitude corrected to z = 2 IShen et al ,l{20 1 1 ). 

c Observations are in decimal years. Dates in parentheses indicate a weighted average of the years preceding it. 

^ Represents the maximum, shared wavelength coverage (rest frame) between two spectra in each source. 

e Indicates whether absorption line variability was detected with the 8 a or 6 <j criterion [yes (Y), no (N), or tentative (T)]. See §3.2 for details. 
f Radio loudness (i.e. R = /6cm//250o) taken from IShen et all (|201 il l ■ 


times ranged between 2 and 3.5 h for each source in the 
sample. 

Observations taken at MDM observatory were reduced 
using standard routines within the image reduction and 
analysis facility (iraf jf|. Processing steps included subtract¬ 
ing off bias structure, flat fielding the data to remove pixel- 
to-pixel variations in sensitivity, extraction of spectra and 
their associated error arrayfl, binning the dispersion axis 
into wavelength units, and applying an extinction correc¬ 
tion and flux calibration to the spectra using a standard 
star from the same night. While the relative flux calibra¬ 
tion of our spectra shows good reproducibility, the absolute 
flux calibration is less reliable in some cases due to variable 
observing conditions. Multiple exposures of a given object 
were combined using an average weighted by the inverse of 
the normalized variance of each exposure. 

Spectra from HET were ac quired using the low- 
resolution spectrograph (LRS; iHill et aD 1 19981 ). Observa¬ 
tions were carried out at an average airmass of ~1.2, with 
exposure times ranging from 660-1920 s. The data were pro¬ 
cessed using standard reduction procedures and were not 
flux calibrated. 

To facilitate comparisons of narrow lines, spectra of a 
given source were smoothed to match the lowest resolution 


Table 2. Spectral parameters. 


Instrument 

A-coverage (A) 

Pixel width (km s *) 

Resolution (km s x ) 

DR7 

3800-9200 

70 

150 

DR12 

3600-10400 

70 

140 

OSMOS 

3100-6800 

40 

190 

CCDS 


80 

200 

LRS 

4250-7400 

50 

230 


Pixel widths and FWHM resolutions represent averages over the 
measured wavelength intervals. Also note that CCDS observa¬ 
tions spanned 1600 A in a single exposure and the coverage win¬ 
dow varied between sources. 

observation via convolution with a gaussian of appropriate 
width. Table 2 shows spectral resolutions representative of 
each instrument used to collect data for our sample. The 
wavelength calibrations were checked using the OI A5577 
atmospheric line and spectra from our entire dataset were 
adjusted as necessary to achieve agreement. We also resam¬ 
pled all spectra to a standard width of 70 km s -1 per pixel 
for our analysis due to varying dispersions between instru¬ 
ments (see Table 2). 


5 IRAF is distributed by the National Optical Astronomy Observa¬ 
tory, which is operated by the Association of Universities for Re¬ 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

6 Extracted error spectra account for CCD read noise and photon 
counting statistics. 


3 ANALYSIS 

In order to study absorption line changes across two or more 
observations in a given source, our general procedure in¬ 
volves two steps: 
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(1) Scale non-DR7 spectra to match the DR7 spectrum 
at non-variable wavelengths using the ratios between the 
observations. 

(2) Quantify absorption line variability by taking flux dif¬ 
ferences between every pair of available spectra and compar¬ 
ing those differences with propagated errors. 

3.1 Scaling spectra 

Spectra were scaled by first computing the ratio between 
DR7 and non-DR7 observations and inspecting the resulting 
ratio to determine which wavelength intervals exhibited no 
‘features’. Features are defined here as abrupt (i.e. occurring 
over intervals < 100 A) changes in the gradient of the ratio, 
and are distinguishable from both the noise and global slope 
changes over a large portion of the wavelength coverage (> 
1000 A). Features in the ratio indicate possible absorption 
line variability and we therefore omitted these wavelength 
intervals when scaling the spectra. 

Changes in the gradient of the ratio occurring over 
> 1000 A were considered featureless wavelength intervals 
and used to scale the spectra. They could be due to quasar 
continuum variations, systematic causes (i.e. flux calibration 
issues or loss of light due to atmospheric dispersion), or per¬ 
haps concurrent variability of overlapping BALs that may 
be confused with the two previous points. This last issue is 
not a problem for our sample for two reasons: 

(i) In sources with overlapping BALs there are also wave¬ 
length regions free of absorption and a feature in the ratio 
would therefore be detectable at the onset of the BALs. 

(ii) Spectra that suffer from absorption over the entire 
wavelength coverage exhibit BALs with structure and well- 
defined edges that would show up as features in the ratio. 

After identifying features in the ratio, fits were per¬ 
formed to the featureless intervals collectively over the entire 
wavelength coverage to produce a scaling function. Func¬ 
tions used in fitting ratios that included MDM and DR12 
epochs ranged between constant values, first ordeiQ linear 
splines, and low order cubic splines; HET data called for 
cubic spline functions up to sixth order due to lack of flux 
calibration. 

When applicable, only featureless intervals not associ¬ 
ated with absorption were chosen to serve as constraints 
for the fit to interpolate or extrapolate over absorption 
lines. There were exceptions where nearly the entire wave¬ 
length coverage was contaminated by absorption (J084044, 
J112220, JI12526, J115436) or the lack of flux calibrated 
spectra from HET required fitting with a non-linear func¬ 
tion (J084044, J121442). If, in these situations, wavelength 
intervals containing absorption were featureless as described 
above, we included these regions when fitting in order to 
scale the spectra in a reasonable way. Two sources exhib¬ 
ited features across nearly the entire wavelength coverage 
and required a separate approach for analyzing the variabil¬ 
ity (J120337, J123103). 

The final step involved multiplying the non-DR7 spec¬ 
tra by the functions produced from fitting the ratios, which 
resulted in scaled non-DR7 spectra that coincided with the 

7 The order specifies the number of spline pieces used for the fit. 


DR7 spectrum at wavelength intervals that were featureless 
in the ratios. 

3.2 Quantifying absorption line variability 

After scaling all available spectra, we adopted a quantitative 
criterion to gauge whether variability was statistically signif¬ 
icant when comparing two spectra from a given object. The 
significance of variability was calculated by taking the differ¬ 
ence in flux between two observations, binning the flux dif¬ 
ference spectrum by 280 km s -1 , and comparing each binned 
pixel with the propagated error a in that pixel. The 280 km 
s' 1 bin width (corresponding to 4 pixels sampled to 70 km 
s -1 ) is greater than the lowest resolution in our dataset (see 
Table 2) and is meant to detect variations across a range 
of pixels instead of probing pixel-to-pixel fluctuations that 
may not reflect true variability. 

Our quantitative criterion for significant variability was 
that the flux difference in a given binned pixel had to be 
at least eight times the propagated error value (8er), either 
positive or negative. The propagated error (a) only accounts 
for CCD noise and random photon statistics. The 8cr thresh¬ 
old accounts for the systematic uncertainties associated with 
the scaling process and correcting for resolution differences 
between spectra, and was determined based upon tests em¬ 
ploying alternative fits to ratios along with convolving spec¬ 
tra with a gaussian of varying width. 

Utilizing the conservative criterion outlined above as 
the sole indicator of variability is questionable given the dif¬ 
ficulty in quantifying the systematic uncertainties, particu¬ 
larly in scaling spectra, due to structural differences between 
the ratios. In light of this we established a threshold for ‘ten¬ 
tative’ variability to be six times the propagated error (6cr) 
for each binned pixel and proceeded to examine each source 
individually to place it into one of three categories: Non¬ 
variable (N), tentative (T), or variable (Y). Each source’s 
category is listed in Table 1 and shown in Fig. 1. Each plot 
in Fig. 1 represents one of the sources in our sample with two 
spectra from different epochs plotted on top of one another 
for comparison. The pair of spectra from each source with 
the maximum, shared wavelength coverage over the absorp¬ 
tion lines and the largest variability by eye is shown in Fig. 
1 . 

Objects in category N did not meet either the 8 a or 6a 
thresholds in any binned pixel. Comparisons that did meet 
either 8<r or 6a yet further inspection by eye suggested the 
variations to be systematic were also conservatively placed in 
category N (see notes in §A.3). Sources fell in the T category 
if they met the 6cr threshold in intervals associated with ab¬ 
sorption with no apparent systematic issues (e.g. J112220), 
or met the 8a threshold but the variability could not be at¬ 
tributed to absorption conclusively (e.g. J112526, J142703). 
Objects categorized as Y met the 8a criterion with none of 
the caveats listed for the N or T categories. The existence of 
underlyin g emission lines in AGN sp ectra (e.g. broad Fell 
emission. IVestergaard fe Wilkesll200lh can make some vari¬ 
able cases ambiguous, and we have taken this into consider¬ 
ation in our assessments of the variability. 

For the purpose of identifying lines as an aid in inter¬ 
preting absorption line variability, we generated a list of 
ionic transitions that is overlaid on all plots in each fig¬ 
ure. The list was produced by inspecting all spectra in our 
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Figure 1. FeLoBAL spectra comparisons. The top of each panel lists the object name, variability category, dates and color of the two 
epochs shown, the rest frame time-scale between the two epochs, and the velocity shift of the overlaid line template (negative velocity 
indicates the lines are blue-shifted). For the error arrays associated with the plotted spectra pairs, the one with larger amplitude is 
indicated by the light blue curve along the bottom of each large panel. The small panels show the binned difference spectrum (blue) 
along with the 8<r (orange, solid) and 6a (orange, dashed) propagated noise spectra. The order of the dates in each panel correspond to 
the order of the difference operation. 
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Figure 1 - continued 
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Figure 1 - continued 
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Figure 1 - continued 
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sample and identifying transitions and multiplets associated 
with absorption, and is based on an inve ntory of transitions 
in BAL quasars compiled bv lHall et al.1 f 2002bl ). We repre¬ 
sent multiplets consisting of more than three transitions as 
horizontal bars extending from the lowest to highest tran¬ 
sition wavelength. The list is not intended to identify all 
transitions contributing to the absorption in any source; line 
blending from the many transitions in Fell and other ions 
make it difficult to identify these lines i ndivi dually. The UV 
multiplet notation comes from I Moor e dl95Ch . 

3.3 Associated absorption lines 

Two sources in our sample (J030000 and J123103) show 
detectable associated absorption lines (AALs) coming from 
Mg II AA2796, 2803 a nd Mgl A2852, as shown in Fig. 2 (see 
e.g. IShen fe Menardll2012l for recent work and general infor¬ 
mation on AALs). Following the procedure outlined above, 
variability between two observations was checked by scaling 
one spectrum to coincide with the other; for this analysis 
we only scaled using intervals around the Mgn and Mgl 
AALs, and used either constants or first order linear splines 
as scaling functions. 

We determined whether AAL variability was significant 
using similar statistical criteria as in the main analysis. The 
requirement for significant AAL variability was that the flux 
difference in a given binned pixel had to be at least six times 
the propagated error value ( 6 <r), either positive or negative. 
Similarly, tentative variability required the binned flux dif¬ 
ference to be at least 4<r. As discussed in §3.2, the propa¬ 
gated error er does not account for uncertainties in scaling 
the spectra or errors from smoothing spectra to correct for 
instrumental resolution differences. The 6 <r and 4<r thresh¬ 
olds for AAL variability are lower than the 8cr and 6 a criteria 
used in the main analysis because there is less uncertainty 
in scaling spectra in the vicinity of narrow absorption lines 
that only span ~10 A. The 6 <r threshold was estimated by 
performing tests involving multiple scaling functions, differ¬ 
ent intervals to constrain the fit, and smoothing spectra by 
various widths to quantify the systematic errors. The 40- 
threshold allows for tentative cases to be quantified objec¬ 
tively. 


4 RESULTS 

We have detected significant broad absorption line variabil¬ 
ity from our 80 - criterion in 3 out of the 12 FeLoBAL quasars 
in our sample and categorized them as Y sources (J120337, 
J121442, and J123103). 

The measured variability time-scales ranged from ~0.6 
to 7.6 yr. Absorption lines that appear to have var¬ 
ied come from Fen (UV1,2,3), Fen* (UV62, 63,780 Mgn 
AA2796, 2803, Mgl A2852, and Nin* (UV12,13) along with 
other Fen multiplets and ionic species that are difficult to 
identify (see §4.2 for evidence of this). Other low-ionization 
lines that were identified but were not associated with re¬ 
gions that met the 8 <r criterion included Fell* (Opt 6 ), Fem* 

8 An * after an ion signifies that the transition originated from 
an excited state. 




2775 2785 2795 2805 2815 2825 2840 2850 2860 2870 

Rest Wavelength (A) Rest Wavelength (A) 


Figure 2. AAL variability in J030000 and J123103. Each source 
shows two panels, each depicting the Mg II AA2796, 2803 and Mgl 
A2852 AALs, respectively. The only distinction compared to Fig. 
1 is that the smaller panels show the 6<r (orange, solid) and 4<r 
(orange, dotted) error spectra plotted with the binned flux differ¬ 
ence. 


(UV34,48), and Aim AA1854,1862. One source in our sam¬ 
ple (J031856) enabled us to examine BALs associated with 
Civ AA1548,1550, SiivAA1393,1402 and Alii A1671. 

For the three Y sources we see evidence for absorption 
line variability across entire BAL troughs, portions of BAL 
troughs, and in narrower absorption lines that do not qualify 
as BALs. We calculated the average velocity width of vari¬ 
able regions meeting our 8 <r criterion within the Mgn AA 
2796, 2803 BAL to be ~520 km s -1 ; this average included 
all regions within the Mgn BAL, whether they overlapped 
with the adjacent Mgn broad emission line (BEL) or Fell* 
(UV62,63) absorption. 

We also looked at how Mgn and Fell varied in rela¬ 
tion to one another. In all three of our Y objects we ob¬ 
served Mg II and Fell BALs varying together, either in¬ 
creasing or decreasing in amplitude. J120337 exhibited de¬ 
creases in absorption line strength over ~3.1 yr across two 
observations; J123103 showed similar behavior over ~2.9 yr. 
J121442 showed decreases in absorption strength followed by 
increases in absorption strength over six observations span¬ 
ning ~7.6 yr. The Y sources are discussed individually be¬ 
low. 

Table_l also shows the radio loudness parameter R 
(from IShen et al.l 1201 ll . defined as the ratio of fluxes at 6 
cm and 2500 A) for all available sources, and we observe 
no notable connection between radio loudness and degree of 
absorption line variability. 
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Rest Wavelength (A) 


Figure 3. Ratio between DR7 and MDM spectra for J120337, revealing evidence of variability in both continuum shape and absorption 
line strength. The blue (dotted-dashed) and red (dashed) lines are the two functions that were used to scale the MDM epoch to coincide 
with DR7 (see §4.1), with the red curve optimized to fit longward of ~2550 A. 


4.1 J120337.92+153006.7 

Fig. 3 shows the ratio between the DR7 and MDM epochs in 
J120337. Inspection of Fig. 3 reveals indication of variability 
over nearly the entire wavelength coverage from ~1940-3000 
A. The global change is not a flux calibration issue because 
this source was observed on the same night as other objects 
that do not show any discrepancies in flux of this magnitude 
when compared to DR7. Slit losses from atmospheric disper¬ 
sion are also likely unimportant given that this source was 
observed at very low airmass (i.e. < 1.1). Fluctuations due 
to quasar continuum emission are also unexpected as they 
typically occur over a much larger wavelength interval. Fig. 
1 shows the resulting variability when the blue dot-dashed 
line in Fig. 3 is used to scale the MDM spectrum; variable 
regions meeting the 8er criterion range from ^2090-2850 A. 

It is plausible that the observed variations are due to 
highly blended absorption lines. Fig. 3 also shows localized, 
abrupt changes which coincide with and mimic the struc¬ 
ture of the absorption lines in the DR7 spectrum in Fig. 1; 
the most notable changes in Fig. 3 occur at ~2720, 2540, 
and 2010 A, corresponding to the Mgn AA2796, 2803, Fen 
(UV1), and Fern* (UV48) absorption lines, respectively. 
Lines that appear in the ratio between two spectra and ex¬ 
hibit the same structure as and align with the actual ab¬ 
sorption lines imply that the absorption lines varied over 
the measured time-scales. However, the possibility remains 
that the observed variability from ^1940-3000 A might also 
be due to peculiar continuum changes. 

In order to conservatively analyze variability due to the 
Mg n AA2796, 2803 BAL, the MDM epoch was rescaled to 
subtract off the large-scale deviation in the ratio spectrum. 
The red dashed line in Fig. 3 shows this scale function as it 
was fit to the ratio spectrum and optimized for wavelengths 
longward of ~2550 A; Fig. 4 shows the resulting rescaled 
MDM spectrum plotted against the DR7 epoch. Using the 
rescaled version of the MDM observation, we detect absorp¬ 
tion line variability meeting our 8<r criterion in three in¬ 
tervals between ^2690-2710 A, with the largest variation 
reaching 9.1<r. The average width of the variable intervals 
was ~490 km s -1 . The Mgn BAL absorbs from approxi¬ 
mately 2670 to 2720 A and likely shortward of 2670 A; the 



Rest Wavelength (A) 


Figure 4. Comparison between DR7 and rescaled MDM spectra 
for J120337 after removing a model of the large-scale variability, 
revealing changes in the Mg II BAL. See Fig. 1 for description of 
the panels. 


lower limit is uncertain due to the onset of possible absorp¬ 
tion from Fen* (UV62,63). 

Although it would be useful to subtract out the large- 
scale deviations and study the residual variability < 2550 A, 
the lack of regions to constrain the fit at the blue end would 
produce large uncertainties. We conclude that the detected 
variability of ,1120337 derives from a unknown change from 
either the quasar continuum or absorption lines along with 
weakening of the Mgn AA2796, 2803 BAL over a ~3.1 yr 
time-scale. 


4.2 J121442.30+280329.0 

Observations for J121442 span time-scales from ~15 d to 
~7.6 yr between any two spectra. We did not detect signifi¬ 
cant variability when comparing any two HET epochs from 
2000, therefore we created a combined HET spectrum via a 
weighted average. 

The shortest time-scale of absorption line variability 
meeting the 8<r criterion was ^0.6 yr between 2010 and 2011, 
and the most significant absorption line changes, at ~20 ct 
significance, occurred over the ~5.9 yr period between 2000 
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Rest Wavelength (A) 

Figure 5. Ratio between the scaled 2010 MDM and DR7 spectra for J121442, revealing differing amplitudes of variability from line to 
line. 


and 2010. The observed changes can be attributed to the ab¬ 
sorption lines because inspection of the ratios shows features 
that mimic the series of lines in the actual spectrum. These 
deviations are localized relative to the wavelength coverage 
and therefore cannot be explained by any intrinsic, contin¬ 
uum variability or systematic issues. 


From 2006 to 2010 we see significant decreases in ab¬ 
sorption line depths from nearly all of the transitions across 
the range 2320-2970 A (see Fig. 1). This pattern met our 
8a criterion for transitions involving Fell (UV1,2,3), Fell* 
(UV62,63,78), Mgi A2852, and the Mgn AA2796, 2803 BAL, 
although it is uncertain whether the last region from ~2780- 
2790 A is due to the Mgll BAL or BEL. We can extract ad¬ 
ditional information by taking the ratio between the scaled 
2010 MDM spectrum and the 2006 DR7 epoch in order to 
see how the strength of the variability varies from line to line; 
Fig. 5 shows the resulting ratio. The small regions around 
2625 and 3000 A are approximately unity and thus serve 
as reference points to where there were no changes between 
the two epochs. It is interesting that the amplitudes of line 
ratios within a given Fe II multiplet are not all the same. If 
we assume that the only contribution to the absorption in 
these regions is due to the identified Fell transitions in an 
optically thin regime then we should expect the line ampli¬ 
tudes within a given multiplet to be the same, independent 
of the origin of the variability. The fact that they are not 
implies that there may be extra contributions to the absorp¬ 
tion lines due to either higher excitation state transitions in 
Fe II or other ionic species, or that some Fe II transitions are 
becoming saturated over our measured time-scales. 


Over the next ~0.6 yr from 2010 to 2011 there are in¬ 
tervals of variability ranging from ~2610-2960 A, which in¬ 
clude Fen (UV1), Fen* (UV62,63,78), the Mgn AA2796, 
2803 BAL-BEL, and Mgi A2852. The absorption line depth 
is increasing for these intervals, opposite to what happened 
in the 2.3 yr before. There is no statistically significant vari¬ 
ability over the ~0.5 yr from 2011 to 2012 or the ~0.6 yr 
from 2012 to 2013, although there is evidence of the con¬ 
tinuing trend of increasing absorption line depth over these 


time-scalefl . When comparing the 2010 and 2013 epochs, 
we detect absorption line variability with the 8<r criterion in 
wavelength intervals ranging from ~2390-2970 A. 

We conclude that across the ~7.6 yr baseline for 
J121442 the absorption line variability is complex. The ab¬ 
sorption lines associated with Fell and Mgi A2852 exhibit 
large, significant decreases in strength over 2.3 yr (see Fig. 
1), followed by increases in depth over the next 1.8 yr (see 
Fig. A). 


4.3 J123103.71+392903.7 

We have a time baseline spanning ~0.5-3.9 yr for J123103 
with four observations for comparisons. Fig. B (see addi¬ 
tional online material) shows comparisons between the DR7 
epoch and the other three available spectra. No absorption 
line variability was detected with 8a confidence between any 
two spectra from 2011-2013. The comparison between the 
DR7 and DR12 epochs in Fig. 1 show that variations occur 
from ^1900-3200 A, but do not exist > 3200 A. The coinci¬ 
dence between the identified ionic species with the variabil¬ 
ity in Fig. 1 is evidence that the observed changes are due to 
the absorption lines and not intrinsic in nature or due to sys¬ 
tematic errors. We note that, in addition to the lines in the 
template, other multiplets of Fe II and other ions are likely 
present to contribute to the absorption from ~1900-3200 A. 

Over the ~2.9 yr between the DR7 and DR12 epochs, 
we see decreases in absorption line strength using our 8<r 
criterion between ~2130 and 2690 A, a region where the 
Fen (UV1,2,3), Fen* (UV62,63), and Nin* (UV12,13) mul¬ 
tiplets occur (see Fig. 1). This source varied with a maxi¬ 
mum significance level of 11a. Similar decreases in absorp¬ 
tion line strength are observed near Fern* (UV48) and Fen* 
(UV78,Opt6), although these do not meet our 8a threshold. 
The BAL extending from ~2650 to 2775 A could consist of 
absorption from both Mgll and Fen* (UV62,63), and we see 
variability only in the high-velocity portion of the BAL. 

We conclude that there is no significant absorption line 


9 This trend is shown visually in Fig. A (see additional online 
material). 
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variability from 2011 to 2013, a time-scale of ~1.0 yr, how¬ 
ever J123103 does exhibit absorption line weakening from 
2005 to 2011 at the 8<r level. 


4.4 Tentative sources 

Sources that meet the 6a threshold for tentative variabil¬ 
ity over wavelength intervals associated with absorption 
(J104459, J112220, J112526) may represent real changes 
that do not meet our 8 a threshold. In order to further as¬ 
sess the variability in these objects, we conducted tests us¬ 
ing larger bin widths than the 280 km s^ 1 width used in the 
general analysis and inspected comparisons shown in Fig. 
1; the larger bin widths probe variations over entire BAL 
troughs, providing additional information beyond the gen¬ 
eral analysis which is optimized to study variability over 
small portions of BALs. 

J104459 shows increases in absorption line amplitude 
centered at ^2650 A and 2700 A, although the former wave¬ 
length is the only interval that reaches 6 ct (see Fig. 1). Tests 
involving larger bin widths yield a maximum significance 
level of ~19cr using a bin width of 10500 km s -1 , correspond¬ 
ing to ~70 per cent of the width of the absorption feature at 
~2675 A. Because the variations are localized relative to the 
wavelength coverage and no other intervals show variability, 
the tentative variations in this source are likely real. The 
absorption may correspond to a high velocity component of 
Mg II. 

The degree of variability in J112220 is significant up to 
~28cr, using a bin width of ^10850 km s” 1 centered on 2590 
A (see Fig. 1); J112526 exhibits 8a variations that may be 
due to the Mg II BEL and we therefore could not obtain a re¬ 
liable estimate of the tentative absorption variability. Both 
J112220 and J112526 do not exhibit localized variability but 
show changes scattered across a large portion of the wave¬ 
length coverage. These variations may be due to Fe II BELs 
instead of absorption lines and the reality of the absorption 
variability in these two sources is therefore less likely. 

,1142703 exhibits 8 a variations from ~2480-2510 A, a 
localized region that is not associated with any transition 
from the line template (see Fig. 1). The variability is un¬ 
likely a flux calibration issue since J121442 was observed on 
the same night as J142703 and does not exhibit similar be¬ 
havior. Variability due to Fe II BELs is also unlikely since 
we would expect changes over larger wavelength intervals. 
The observed variation is therefore potentially due to ab¬ 
sorption and might be associated with highly blue-shifted 
Mg II AA2796, 2803. ,1084044 is not included in this analysis 
due to the ambiguity in associating variations with absorp¬ 
tion features. 

We conclude that, of the sources categorized as tenta¬ 
tive, J104459 and J142703 exhibit absorption line variabil¬ 
ity that is likely real and associated with blue-shifted Mg II 
AA2796,2803. 


4.5 Associated absorption lines 

Fig. 2 shows spectra comparisons of the Mgn AA2796, 2803 
and Mg I A2852 AALs for J030000 and J123103, the only two 
sources in our sample that have AALs. J030000 did not show 
variability under our 6<r criterion between any of the three 


DR7 observations; we therefore used the weighted average 
epoch to compare with the MDM observation, a time-scale 
of ~6.0 yr. 

We detect tentative variability in ,1030000 meeting the 
4a criterion in the Mgn A2803 AAL and also at ~2815 
A, which may be associated with red-shifted Mgn or some 
other unidentified ion (see Fig. 2). Both of the Mg Hand Mgl 
AALs increase in amplitude upon inspection by eye, how¬ 
ever the variations occur across ~2 pixels, making it hard to 
consider the changes as real. It is interesting to note, how¬ 
ever, that the Mgn AALs are at equal depths but do not 
reach zero flux, an indication that Mgn is saturated and 
partially covers the continuum source. Using our 6er thresh¬ 
old, we conclude that there is no significant AAL variability 
in J030000. 

We also did not detect AAL variability at 6 a between 
the four observations of J123103, with time-scales up to ^3.9 
yr. We note, however, that the comparison between 2011 
and 2013 (see Fig. 2) shows tentative changes at 4cr in the 
Mgl A2852 line. The amplitude of the Mgl AAL appears 
to increase over the ~1.0 yr time-scale and the variation 
persists over multiple pixels, however we conclude that the 
changes are only marginal. 

The lines we detect in ,1030000 and ,1123103 are likely 
composed of narrower, unresolved features; average mea¬ 
sured FWHM for the Mg II and Mg I AALs are ~150 and 135 
km s' 1 for J030000, and ~270 and 210 km s -1 for J123103, 
respectively. We note that a portion of the AAL profiles in 
J123103 is associated with blue-shifted absorption, while the 
rest appears to be red-shifted. This behavior, combined with 
the potential red-shifted Mgn AAL system in J030000, indi¬ 
cates that the AAL gas is complex with possible outflow and 
inflow components. We conclude that there is no AAL vari¬ 
ability meeting our 6er criterion in either of the two sources 
analyzed, however future observations may provide insight 
into the nature of the possible variations. 

5 DISCUSSION 

5.1 Interpreting the variability 

Placing constraints on outflows requires an understanding 
of the physical origin of the absorption line changes. Two 
scenarios that can potentially explain the variability are that 
either the absorbing gas undergoes an ionization change or 
a change in covering fraction occurs along our LOS due to 
transverse motions of gas across the emission region of the 
AGhEJ Below we discuss the origin of the variations on 
the basis of photometric and optical depth arguments, and 
provide evidence that broadly supports change in covering 
fraction as the preferred scenario to explain absorption line 
variability in our sample. 


10 A variant of the second scenario is that a portion of the con¬ 
tinuum source outside the absorber sightline undergoes a local 
change in emissivity, leading to fluctuations in an unabsorbed 
continuum component. 
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5.1.1 Photometric arguments 

The photometric behavior of our sources may provide in¬ 
formation regarding the origin of the absorption line vari¬ 
ations. We utilize synthesized CSS V-band photometry to 
look for significant flux variations over our measured time- 
scales. Fig. 6 shows photometry of the three Y sources along 
with spectral observation dates from our sample (red dot¬ 
ted lines fH. There is a small decrease (~0.16) in V-band 
magnitude between the two observations in J120337; a de¬ 
crease of ~0.18 mag is observed between the first two epochs 
in J121442, which is followed by a smaller increase from 
^2010-2012. Photometry provided for J123103 is less reli¬ 
able due to large variations that likely originate from sys¬ 
tematic errors (the source is close to a bright star). 

Comparisons between our spectroscopic findings and 
measurements in Fig. 6 reveal a noteworthy result in that 
absorption line strength correlates inversely with estimated 
V-band flux. This behavior is seen for each of the epochs 
in J120337 and J121442 and possibly in J123103. A quali¬ 
tative pattern of this type is expected in scenarios involving 
gas absorbing in the V-band and crossing the LOS while the 
background light source remains constant. Alternatively, the 
V-band flux primarily traces changes in the source contin¬ 
uum luminosity, with the absorber ionization state and ion 
column density responding to the incident radiation field to 
produce the observed trend. Since we do not have absolute 
photometry for most spectra, it is ambiguous whether the 
background light source changed in our objects. However, a 
calculation of the fractional change in synthesized V-band 
flux for our scaled spectra can be compared with the CSS 
measurements to constrain the intrinsic continuum variabil¬ 
ity. 

The spectra shown in Fig. 1 would thus predict for 
J120337 a ~40 per cent increase in synthesized V-band flux 
between 2005 and 2012; the CSS measurements show a much 
smaller increase of ~13 per cent. The difference can be un¬ 
derstood if the weakening of absorption is accompanied by 
a drop in continuum luminosity, however inspection of Fig. 
1 suggests an alternative explanation of a change in con¬ 
tinuum shape or possibly systematic errors. For J121442, 
our wavelength coverage does not allow us to synthesize 
a V-band flux, however using spectra from 2006-2010 we 
compute a ~23 per cent increase in flux over 3800-5150 A; 
this is slightly larger than the ~15 per cent change derived 
using CSS fluxes near those dates. If the modulations are 
only due to absorption and variable coverage, we would ex¬ 
pect higher amplitudes at shorter wavelengths in this source 
due to the larger number of transitions, consistent with this 
comparison. Our third source, J123103, yields a ~20 per 
cent increase in V-band flux between 2005 and 2011, calcu¬ 
lated using spectra from our sample. CSS photometry for 
J123103 has systematic uncertainties and therefore a com¬ 
parison could not be tested quantitatively, however there is 
also evidence of a flux increase from 2005-2011. 

In summary, with the current data we cannot distin¬ 
guish on the basis of photometric arguments whether vari- 


11 Photometric changes in the other sources in our sample were 
found to be < 0.1 mag when comparing CSS measurements at 
the epochs for our spectroscopic observations. 
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Figure 6. Estimated CSS F-band photometry for the Y sources, 
depicting variations between spectral observations (shown as red 
dotted lines). Each upper panel shows all photometric measure¬ 
ments with their errors, and each lower panel shows average mea¬ 
surements over one year increments. Blue symbols above the red 
dotted lines indicate the directional change in absorption line am¬ 
plitude relative to the previous spectral observation; absorption 
lines either decreased (-), increased (+), or did not meet our 8cr 
criterion (0). 


ability is more likely due to a change in coverage or an ion¬ 
ization change. 


5.1.2 Optical depth arguments 

Variability in absorption lines that are known to be optically 
thick, if observed, provides strong evidence for changes in 
coverage (e.g. Capellupo et al. 2014a). The Mgn feature can 
potentially be used as such a diagnostic, and information 
on its optical depth can be derived from considering the 
combined behavior of Mg II and Fe II absorption. 

For that purpose, we calculated a different optical 
depth i~Fe for each Fell multiplet present in our sample 
(UVl,2,3,62,63,78,Opt6) and use those values with the op¬ 
tical depth TMg, calculated using the Mgn AA2796, 2803 
doublet, to determine a set of ratios TMg/rFe that depend 
on the respective ionic column densities, transition oscil¬ 
lator strengths, and transition wavelengths. For each ra¬ 
tio we adopted ave rage transition wave lengt hs and oscilla¬ 
tor st rengths from IVerner et al.l |l996) and iFuhr fe Wiesel 

mi). 

For our calculations, we assumed the Mg II to Fe II col¬ 
umn density ratio to be equal to the solar elemental abun- 
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dance ratio (taken from lAsplund et al .1 l20()il ). We did not 
factor in the unknown level populations of Fe II and Mg II or 
the distribution of iron and magnesium over their ionization 
stages, but adopted the solar abundance ratio throughout. 
Incorporating details about the physical conditions of the 
gas, however, only favors larger TMg/rFe than the ones cal¬ 
culated here due to the large number of low-lying Fe II energy 
levels and closely spaced ionization potentials in iron relative 
to magnesium. Since we are only concerned with establish¬ 
ing very crude lower limits on TMg/TFe for the purpose of 
understanding the variability, the above assumptions do not 
affect our interpretations. 

We calculate lower limits on TMg/rFe to be, on aver¬ 
age, ~6.3 using Fen (UV1,2,3) multiplets, ~4.1 using Fell* 
(UV62,63), ~15 using Fen* (UV78), and ~50 for Fen* 
(Opt6). Observing Fen in absorption at a given optical 
depth therefore implies that the Mgn optical depths can 
range from being at least ~4.1-50 times larger than Fell, 
depending on which Fell multiplets exhibit absorption. We 
use this information below to assess whether a change in cov¬ 
ering fraction due to transverse motions of gas is the likely 
scenario occurring in the three Y sources. 

J120337 appears to exhibit absorption from Fell* 
(UV62,63) at low optical depths but no detectable absorp¬ 
tion from Fen* (UV78,Opt6), and Mgn is therefore less 
likely to be saturated (see Fig. 1). We concluded in §4.1 
that the observed variability in this source is due to both 
an unknown change across the entire wavelength coverage 
and a weakening in the Mgn AA2796,2803 BAL. Localized 
changes in the ratio spectrum are detectable near transi¬ 
tions from Fell (UV1) and Fern* (UV48), implying that 
these lines may have varied in the same manner as Mgn 
(see Fig. 3). However, the poor signal-to-noise ratio towards 
the blue end and the fact that we cannot fit these regions 
with any degree of certainty makes it inconclusive if this 
behavior is occurring in a significant way. Ions of different 
ionization state varying in concert supports the scenario of 
a change in covering fraction due to transverse motions of 
gas, however the fact that Fein has a ground state ioniza¬ 
tion potenti al not too dis similar from Fen (~31 and 16 eV 
respectively, iDrain IB makes it difficult to rule out an 
ionization change in this source. 

J121442 shows detectable absorption coming from Fell 
(UV1,2,3), Fen* (UV62,63,78), and Fen* (Opt6); this im¬ 
plies that Mg II is very likely saturated in this source (see Fig. 
1). Significant variability meeting the 8<r criterion occurred 
in regions associated with Fen (UV1), Fell* (UV62,63,78), 
and the Mgn AA2796, 2803 BAL (see Fig. 1), although 
the variable interval near Mgn may be due to either the 
BAL or BEL. The lack of variability in the deep portion 
of the Mg II BAL might be due to differing covering frac¬ 
tions between the Mgn and Fell gas in an inhomogeneous 
absorber; large TMg/rFe ratios along with closely spaced ion¬ 
ization potentials in iron relative to magnesi um s u pport 
a hig her Mgn cove ring fraction than Fe n (see iLucv et al.l 
l2014h . Futhermore. Ide Kool et all (l2002al l used a template¬ 
fitting approach to model the spectrum of J121442 and con¬ 
strained the Fen column density to be ^io 16 ' 95±01 ° cm -2 . 
They concluded that the high column densities imply that 
Fell is saturated and thus partially covers the continuum 
source; Fig. 1 shows evidence of this as the regions where the 
Fell (UV1,2,3) multiplets exist are approximately flat and 


do not reach zero flux. The likely interpretation to explain 
the variability in J121442 is therefore a change in covering 
fraction from transverse motions of the outflowing gas. 

Mg II is also likely saturated in ,1123103, as this source 
exhibits absorption coming from Fell (UV1,2,3) and Fell* 
(UV62,63,78), with the possible presence of Fell* (Opt6) at 
very small optical depths. Significant decreases in absorption 
line strength meeting the 8<r criterion are observed in the 
regions with transitions from the Fell (UV1,2,3) and Fell* 
(UV62,63) multiplets; no variability is detected in the lower 
velocity portion of the Fen* (UV62,63)-Mg II AA2796, 2803 
BAL trough (see Fig. 1). The same changes are observed 
near Fern* (UV48) and Fen* (UV78,Opt6), although these 
regions did not meet our 8<r criterion. 

An ionization change within the gas that ionizes Fe II to 
Fern would predict weakening in Fell absorption lines and 
strengthening in the Fern, which is not observed. It is possi¬ 
ble that both ions are decreasing in favor of higher ionization 
stages. It is noteworthy, however, that only a portion of the 
Fern trough weakens. Detecting variability in portions of 
BALs can be explained by a change of covering fractions due 
to the observed ionic species moving out of the line of sight 
only at the higher velocities; evidence of saturation in the 
Fell (UV1,2,3) and Nill (UV12,13) lines also lends support 
to transverse motions of gas in J123103, as the wavelength 
intervals across these multiplets are approximately flat and 
do not reach zero flux (see Fig. 1). If Mgn is saturated, how¬ 
ever, the lack of variability in the lower velocity portion of 
the BAL around 2750 A is also consistent with an ionization 
change scenario. It is possible that this source is undergoing 
both gas motion across the LOS and ionization changes; Fig. 
1 also shows weakening in the Fen* (UV78) line across its 
entire absorbing range. 

In light of the above evidence for each of our Y sources, 
the observed changes are readily interpreted as a change in 
covering fraction, and this explanation is preferred in partic¬ 
ular for J121442 and J123103. We therefore consider an out¬ 
flow crossing the LOS in order to place physical constraints 
on the outflows in our sample. 

5.2 Constraining FeLoBAL outflows 

Evidence favoring transverse motions of gas across the LOS 
to explain the observed variability allows us to constrain 
the distance of the absorbing ions from the SMBH using the 
measured variability time-scale At, an estimate of the size 
of the continuum region at a particular wavelength, DA peak , 
and the approximate fraction A A of the continuum region 
that the outflow traversed over the measured time-scale. 
Following ICapellupo et ahl (1201 lh . we quantify A A by mea¬ 
suring the average flux of a variable region for two epochs 
of interest (/i and / 2 ) along with the average continuum 
flux f c across the same interval; we then calculate the ratio 
(/i — fi)/fc to serve as the estimate of A A. Due to the many 
blended BALs present across our sample, determining where 
to place the continuum was difficult and highly uncertain. 
We therefore conducted multiple tests, placing the contin¬ 
uum at extreme locations that we believe should encompass 
the range of possible continuum values, and quantified aver¬ 
age estimates on A A and its associated error for each source. 
Our estimates for A A were 0.21 ± 0.10, 0.19 ± 0.11, and 
0.34 ±0.18 for J120337, J121442, and J123103, respectively. 
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For an estimate of LU peak , we utilize the result from 
Chen et al. (in preparation) who assume a standard thin ac¬ 
cretion disk that emits as a blackbody and use Wien’s dis¬ 
placement law to translate the temperature at a given radius 
to the wavelength corresponding to maximum blackbody 
emission to estimate the radius of the disk at that wave¬ 
length. Their result adopts a modified constant in the radius 
function derived using the above assumptions; the new fac- 
tor is meant to comprom ise between this standard analysis, 
iBlackburne et all (120111 1 who us ed micro-lense d quasars to 
estimate accretion disk sizes, and lGaskelll (120071 ') who argued 
that the temperature profiles of real AGN disks are flatter 
than the standard predictions. 

The calculations of D A peak require estimates of the 
Eddington ratio Tboi/Tedd, SMBH mass M ,, radiative ef¬ 
ficiency 77 , and maximum blackbody emission wavelength 
Apeak- We took values for Lboi/L e dd and M . from lShen et al.l 
teOllIf to be ~0.06 and 7 x 10 9 M Q for J120337, ~0.07 and 
3 x 10 9 Mq for J121442, and ~0.03 and 4 x 10 9 M Q for 
■112310(1^1. We did not use their quoted statistical errors on 
Tbol and M, for our calculati ons; we inste ad adopt uncer¬ 
tainties based on estimates bv IShenl (120131 ') to be ~0.3 and 
0.5 dex for Lbol and M ,, respectively, reflecting syste matic 
error s associated with single epoch measurements (seejShenj 
120131 for a review of quasar black hole mass measurement 
methods and discussion of these uncertainties F 9 !. We adopt 
0.10 as our estimate of 77 , and calculate an average value of 
Apeak to be ~2700 A based upon the location of the variable 
regions in the three Y sources. 

Using the estimates for the respective continuum source 
sizes, we consider the extreme case of a circular outflow of 
infinite radius crossing a squar e continuum (the so-c alled 
‘knife-edge’ model taken from ICapellu po et all 120131, see 
their f ig. 14b). We follow the calculation of Capellupo et al.l 
(|2013i ) with this model to determine the crossing speed to 
be Derosa = AA x L>A peak /Af. We can then tie the crossing 
speed to a distance assuming that the outflow is in keplerian 
rotation around the SM BH (i.e. r = GM, / v^ rnsK ). The knife- 
edge model presented by lCapellupo et ~akl (120131 1 is meant to 
be a limiting case of all geometrical situations involving an 
outflow traversing our LOS to the continuum-emitting ac¬ 
cretion disk, thereby establishing a lower limit on Dcross and 
upper limit on the distance r from the SMBH. 

Using the above parameters and their errors, we 
estimate diameters of the continuum emitting regions 
for J120337, J121442, and J123103 to be ~0.031±£;^, 
0.018 !q;qJ°, and 0 . 0161 q;qq 9 pc, respectively; the errors were 
propagated throug h using the estimated uncertainties on 
Lbol and M. from IShenl (|2013l ). Using our respective esti¬ 
mates on A A and their errors, we constrain lower limits on 
the crossing speeds to be Derosa > 550, 990, and 540 km 
s -1 for J120337, J121442, and J123103, respectively. These 
yielded upper limits on the outflow distance from the SMBH 


12 Single epoch estimates of M, for our Y sources were based on 
either the Mgn (J120337, J123103) or H /3 (J121442) BEL. 

13 Note in particular that FeLoBAL spectra may typically be 
more reddened than the rest of the quasar population (e.g. 
ISpravberrv_fc_Folt 2 ^ 1992|b therefore bolometric corrections de¬ 
rived bv IShen et al. l201ll ') may not be representative of these 
sources and can lead to systematic errors in Lbol and M.. 


to be r < 26, 2.5, and 23 pc for our three Y sources, respec¬ 
tively. 

The upper limits on r do not account for errors on 77 
and Apeak- While A pea k can only fluctuate on the order of 
tens of angstroms over our variable regions, 77 can reach up 
to 0.42 (for a maximally spinning Kerr black hole in co¬ 
rotation with the accretion disk). Our adopted 77 ~ 0.1 es¬ 
timate is potentially uncertain, given that recent work has 
estimated high spins and radiative effi ciencies on massive 
black holes in quasars from z ~ 1.5-3.5 (lTrakhtenbrot|[2014l : 
ICapellupo et al.ll2015T) . Since we do not have information 
about the spins of the SMBHs in our sample, we consider 
the extreme case by adopting 77 to be 0.42 and A pea k to be 
the shortest wavelength in each variable region used; these 
extreme parameters yield the maximum possible upper limit 
on r for each of our variable sources. When also accounting 
for errors from the other parameters, we thus derive conser¬ 
vative upper limits to be r < 69, 7, and 60 pc for J120337, 
J121442, and J123103, respectively. 


5.3 Comparisons to previous work 


Our measured variability time-scales of ~0.6-7.0 yr and our 
limits on outflow distance for the Y sources are broadly 
consiste nt wit h prev i ous v ariability studies of FeLoBAL 
quasars. IVivek et al.l d 2012 h saw variable time-scales down 
to ~1.2 yr in one out of the five FeLoBAL quasars in 
their sample, and attri buted the chang es to Fern and 
Aim fine-structure lines. [Hall et alJ (i201ll ) detected absorp¬ 
tion line variability over 0.6-5 yr in the FeLoBAL quasar 
J1408+3054, observing dramatic changes in Mgll and dis¬ 
appearance of Fell lines at the same velocity; using a dif- 
ferent outf low model than the one adopted in our analysis, 
lHall et all (120111 ) constrained the absorber distance to be 
between 1.7 and 14 pc from the SMBH, in approximate ac¬ 
cord with our results. 

The general consensus from previous work on HiBAL 
quasars is that BAL variability is common over multiple 
year time-scales (see iFiliz Ak et al.l 120131 . ICapellupo et al.l 
l 20 ll and refs, therein). This is consistent with our results. 
Previous studies on C IV and Si IV BAL variability have also 
concluded that variable regions tend to occur in portions 
of BAL troughs dGibson et al.ll20 0Sl: ICap ellupo et ahll2012l . 


120131 : IFiliz Ak et al. 20131 1. Filiz Ak et al.l ( 20131 ') calculated 
the mean velocity width of variable regions across their sam¬ 
ple of Civ and Siiv BALs to be 713.6 km s _1 and 592.8 km 
s -1 , respectively. These are on the same order of magnitude 
as our Mg 11 BAL estima t e of <^520 km s -1 for our three Y 
sources. I Capellupo ef~ai1 d201ll. 120121. l2013lf constrain their 
sample of variable C iv and Si iv BAL outflows to distances 
< 10 pc. 

Previous work that has constrained quasar outflow dis- 
tance s using p hot o ioniza tion modeling yields a range of re¬ 
sults. iDunn et ai] Unci) studied J0318 56, categor i zed as a 
‘bright, reddened’ FeLoBAL quasar bv lHall et al.1 d2002bl l. 
and divided the BALs into approx i mate ly eleven separate 
kinematic components. IPunn et al.l ( 2010 ) also determined 
the electron density of the gas using metastable excited lev¬ 
els from Sill, Fell, and Nill and, combined with an estimate 
of Lbol, constrained the distance of the outflowing gas to 
likely be ~6 kpc from the SMBH. Our lack of variability 
in J031856 with our available data is consistent with the 
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distance estimate of [Dunn et ali J2010l i: photoionized gas 
at larger distances from the SMBH implies lower densities 
(for a given ionization parameter) which result in longer re¬ 
combination time-scales for an ionization change to occur, 
and transverse motions across the line of sight are less likely 
to produce detectable changes from absorbing ions farther 
awa y from the SMBH. 

lArav et all (i2013h used density-sensitive, unblended ab¬ 
sorption lines from O IV in the far UV spectrum of the 
non-BAL quasar HE0238-1904 to estimate an o utflow dis¬ 
tance of 3 kpc from the SMBH. Ide Kool et al.l (|200ll '> an¬ 
alyzed the FeLoBAL quasar J104459 by assuming a single¬ 
phase outflow model and determined the distance of the Fe II 
and Mg I absorb in g reg io n fro m the SMBH to be ~700 pc. 
lEverett. Konigl fe Aravl (12002T ) conducted a follow-up anal¬ 
ysis of this source, instead assuming a shielded, multi-phase 
outflow that consisted of liigh-density clouds embedded in a 
low-density, continuous magnetohydrodynamic wind. This 
reinterpretation to explain the observed nature of the ab¬ 
sorption lines in J104459 yielded a new distance estimate be¬ 
tween the absorbing gas and SMBH to be ~4 pc. These two 
estimates are model dependent and differ by approximately 
two orders of magnitude. Over our available time-scale of 
~4.0 yr we detect tentative absorption line variability with 
our 6 <t criterion for J104459. It is interesting that the tenta¬ 
tive variations are likely from a localized region that could 
be highly blue-shifted Mgn, while the other absorption line 
systems do not vary. Our resu lts are therefore broadly con¬ 
sistent with both Ide Kool et al.l (1200il l and Everett et al. 
( 2002 ) in the sense that this source may have outflowing gas 
at two distinct radii. 

Order of magnitude differences in the distances of out¬ 
flows from the SMBH have direct implications for theoret¬ 
ical models attempting to explain the outflow origin and 
underlying acceleration mechanisms. Quasar outflows have 
been modeled as winds arising in the same regio n as the 
BELs on sub-parsec sca l es (se e e.g., Murrayjetld] 


Progg : i _Stone^^KaIhnar] l2000l . IProea fe KallmaiJ 


1995, 


2004, 


Fukumura et al • 20 ld), while Faucher-Giguere et al. (2012) 


explained FeLoBAL formation in particular to be from the 
quasar wind propagating out to kiloparsec scales and in¬ 
teracting with ISM gas to form clum ps tha t prod uce the 
obser ved absorption (see also iFaucher-Giguere fe Quataertl 

l2012jb _ 

iLucv et alJ (2014) discussed the selection bias associ¬ 
ated with photoionization models that require unblended 
lines from both ground and excited states of a given ion 
to constrain the density of the outflowing gas; they argue 
that density-sensitive lines observed in BAL quasar spectra 
are only useful in lower density gas and outflow distance 
constraints are therefore larger. Looking for objects with 
these characteristics may inherently bias the result and sub¬ 
sequent con clusions may not e xtend to the entire population 
of outflows. ILucv et all (! 2014l l estimated the outflow in the 
FeLoBAL quasar FBQS J1151+3822 to be between 7.2 and 
127 pc from the SMBH using photoionization modeling and 
blended, excited state lines in Fell and Hel, a range that is 
similar to results from variability studies. 

Our upper limit of ~7 pc for_ J121442 is consistent with 
the results of Ide Kool et al.l ri 2002 aQ, who studied a Keck 
HIRES spectrum of this source by fitting a model to the ob¬ 
served data in order to extract column densities, excitation 


temperature, and additional information on the continuum 
and covering factor of the outflow. In co njunction with pho¬ 
toionization models. Ide Kool et al . ([2002a.) were able to con¬ 
strain the distance of the outflow in J121442 to be between 
1 and 30 pc from the SMBH. These distance constraints 
should also broadly apply to other quasars that have similar 
spectra as J121442 (i.e., large Fen column densities and the 
presence of Fell* absorption lines); sources in our sample 
that have similar luminosities and absorption spectra, and 
therefore are appro ximately consis t ent wit h the 1 — 30 pc dis¬ 
tance constraint of lde Kool et al] d2002al l. include J030000, 
J084044, J112220, J112526, J115436, and J123103 (see Fig. 
!)■ 

Other studies have fo und a range of outflow distances to 
exist in the same object. Ide Kool et al.l d 2002 bl l analyzed a 
Keck HIRES s pectru m of J084044, a ‘radio-moderate’ BAL 
quasar dBecker et al.lfl99^h . concluding that there are two 
distinct kinematic outflow components with different den¬ 
sities; the distances of the high and low-density absorbing 
regions from the SMBH were found to be ~1 and 230 pc, 
respectively, based on column density measurements of Fe ill 
and Al ill along with level popu lations obse r ved in low-lying 
states in Nill, Sill, and Fell. iMoe et ail d2009l f analyzed 
the FeLoBAL quasar J0838+2955, constraining the lower 
velocity BAL system to be at 3.3 kpc from the SMBH us¬ 
ing Sin* and Fen* lines with photoionization modeling but 
also detected variability in the higher velocity C IV BAL sys¬ 
tem and_constrained it to be < 0.1 pc from the SMBH. 
iFiliz Ak et aid (|2013i f detected both variable and stable 
CIV troughs over their measured time-scales in two sources 
(J141513.99+365412.2 and J225608.48+010557.7, see their 
§5.1). In light of the above discussion, it seems plausible that 
FeLoBAL outflows can exist at a wide range of distances 
from the SMBH. 


6 SUMMARY 

We have detected absorption line variability using our 80 - 
criterion in 3 out of the 12 FeLoBAL quasars in our sample 
over rest frame time-scales of ~0.6-7.6 yr (J120337, J121442, 
J123103). We have also further quantified and assessed the 
significance of variability in our tentative (T) sources, con¬ 
cluding that two objects likely exhibit real variations in 
the absorbing region (J104459, J142703). We have also an¬ 
alyzed Mg II and Mgl AAL variability in two of our sources 
(J030000, J123103), finding only marginal variability using a 
tentative 4a criterion over measured time-scales up to ~3.9 

yr- 

Wavelength intervals of variability that met our 8 er cri¬ 
terion are associated with transitions coming from ground 
and excited states in Fe n as well as Mg n AA2796, 2803, Mg i 
A2852, and Nill* (UV12,13). Due to our limited wavelength 
coverage and signal-to-noise ratio, we were unable to probe 
variability of Fe ill and higher ionization species such as C IV 
and Si IV. Studying a larger range of ionization states simul¬ 
taneously would be useful in the future for interpreting the 
origin of the observed variability. 

Sources meeting our 8a threshold (J120337, J121442, 
J123103) exhibit absorption line changes readily interpreted 
as transverse motions of gas across the LOS; this scenario 
is preferred in particular for J121442 and J123103 due to 
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evidence of saturation in the Fell, Nill, and Mgn absorp¬ 
tion lines. Other scenarios cannot be entirely ruled out. In 
the crossing absorber scenario, we constrain the outflow dis¬ 
tances from the SMBH to be < 69, 7, and 60 pc for J120337, 
J121442, and J123103, respectively, with conservative as¬ 
sumptions. 

Our results differ by a few orders of magnitude with 
some studies that utilize photoionization models and con¬ 
strain outflows to be at kiloparsec distances from the SMBH. 
Two possibilities are that either two outflow phases of differ¬ 
ent origin are being probed, or we are tracing the same out¬ 
flowing gas in different stages in its evolution. Determining 
the incidence of absorbing outflows as a function of distance 
from the SMBH using large sample studies will be important 
for understanding this possible connection. 
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APPENDIX A: NOTES ON INDIVIDUAL 
OBJECTS 

Al J030000.57+004828.0 

There are three available DR7 observations for J030000; to 
compare these epochs we scaled the two most recent spectra 
to coincide with the earliest observation. We did not detect 
any significant variability with the 8<r criterion between any 
two DR7 epochs so we combined the three spectra together 
via a weighted average to be compared to the MDM observa¬ 
tion; no significant variability was detected. The time-scales 
bet ween any two epoch s ranged from ^35 d to ~6.0 yr. 

IVivek et al.1 (12012!) also conducted a variability study on 
this source, utilizing observations in 2007 and 2008 in addi¬ 
tion to the DR7 observations. They did not detect significant 
BAL variability across their ~4.2 yr coverage, consistent 
with our results, however they report a possible variation 
of the feature between ~2400-2500 A. They see a consis¬ 
tent rise in flux across this interval between DR7 and their 


2008 spectra, however we do not see any change in our com¬ 
parison between DR7 and 2012 (see Fig. 1). It is unclear 
which io nic species cont ribute to the absorption in this re¬ 
gion; IVivek et al.l (l2012l l mention the possibility of a higher 
velocity component of the Fell (UV1) multiplet. The lack of 
variability that we observe between 2000 and 2012 suggests 
the ob ject varied over the time-scale probed bv l Vivek et al.l 
J2012ll then went back to its original state; follow-up obser¬ 
vations may shed light on this behavior. 

A2 J031856.62-060037.7 

We have three observations for J031856 spanning time-scales 
from ~3.4-3.8 yr. The wavelength coverage for studying 
variability was from ~ 1300-2030 A, a range that included 
the resonance transitions Siiv AA1393, 1402, Civ AA1548, 
1550, Al 11 A1670, and Al 111 AA1854, 1862. The Fe 11 and Mg 11 
transitions are found in redder regions of the spectrum out¬ 
side of our observational coverage. We do not detect any sig¬ 
nificant absorption line variability using the 8cr criterion for 
these ions ac ross the obs e rved time-scales, consistent with 
the results of lVivek et all (12012! ) who also study this source 
over ~3 yr between 2001 and 2009 and had additional wave¬ 
length coverage extending out to ~2850 A. 

A3 J084044.41+363327.8 

Our sample contains eight spectra for J084044, spanning 
time-scales from ~10 d to ~5.3 yr. No variability was de¬ 
tected between the HET epochs using the 8cr criterion so 
we combined them into one weighted-averaged spectrum for 
comparisons (weighted averages were also employed for HET 
epochs of J121442 and J142703, due to lack of variability). 

Possible absorption line changes were detected over our 
available time-scales in scattered regions between ~2210- 
2800 A. The region at ^2800 A cannot be conclusively at¬ 
tributed to the Mg 11 AA2796, 2803 BAL due to contamina¬ 
tion with the Mg 11 BEL. In addition, based upon follow-up 
inspection, we conclude that the variability elsewhere is po¬ 
tentially due to systematic uncertainties in matching instru- 
rne ntal resolutio ns._ 

IVivek et al.l (12012! ') conducted a variability analysis on 
J084044 from 2001 to 2009 (~4 yr), finding no significant 
absorption line changes. Our lack of conclusive absorption 
line variability on time-scales of ~0.85-5.3 yr is consistent 
with their result. 

A4 J112526.12+002901.3 

The time-scales between our three available epochs for 
J112526 range from ~0.55-6.4 yr. The only significant vari¬ 
ability was detected between the DR7 and DR12 compar¬ 
ison, a timescale of ~5.8 yr, with variable regions ranging 
from ^2770-2810 A. This could correspond to a weakening 
in the Mgn AA2796, 2803 BAL and/or strengthening of the 
Mg 11 BEL. lHall et al.l (j2013i) came to the same conclusion 
after scaling the DR12 spectrum using a function expressed 
as a constant times a power law, different than the first order 
linear spline function we used to scale the DR12 epoch. 

Inspection of Fig. 1 shows that there are tentative de¬ 
creases in absorption line strength coming from the Fell 
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(UV1,2,3,62,63) multiplets. Because these regions only meet 
the 6 <t criterion and due to the contamination of the Mgn 
BAL with the Mgn BEL, we conclude that J112526 is a 
tentative case and the origin of the variability is ambiguous 
(although see §4.4 for discussion). 

This source is one of the few B AL quasars for which 
a portion of the BAL is red-shifted (lHall et al.ll2013l ). The 
velocity range over which we see Mgn absorption i s from 
approximately —1000 to +1400 km s -1 . lHall et ahl (12013T ) 
also mentioned that the blue-shifted part of the Mg II BAL 
weakened more than the red-shifted part, an observation 
that we also detect by eye. Two likely scenarios to explain 
red-shifted absorption could be either a rotationally domi- 
nate d outflow, radia l infall of gas or a combination of both 
(see [Hall et al.|[2013l for details). Observing a higher degree 
of variability in the blue-shifted part of the Mgn BAL as 
opposed to the red-shifted part might indicate that there 
are multiple outflow-infall systems contributing to the ab¬ 
sorption. 

This paper has been typeset from a TJjjX/ file prepared 

by the author. 










20 S. M. McGraw et al. 



Rest Wavelength (A) 


Figure A. Spectra comparisons of the FeLoBAL quasar J121442, revealing a change in the direction of absorption line variability. The 
top of the figure lists the object name and velocity shift of the overlaid line template (negative velocity indicates the lines are blue-shifted), 
which is shown in the topmost panel. Decimal years next to each comparison indicate the two epochs plotted, with the color of each year 
corresponding to the color of the associated spectrum. Solid horizontal lines represent zero points for each comparison. 
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Figure B. Spectra comparisons of the FeLoBAL quasar J123103. The top of the figure lists the object name and velocity shift of the 
overlaid line template (negative velocity indicates the lines are blue-shifted), which is shown in the topmost panel. Decimal years next 
to each comparison indicate the two epochs plotted, with the color of each year corresponding to the color of the associated spectrum. 
Solid horizontal lines represent zero points for each comparison. 







































































































































